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The competition between intertube hopping processes and density-density interactions is investigated in one- 
dimensional quantum dipolar bosons systems of coupled tubes at zero temperature. Using a phenomeno- 
logical bosonization approach, we show that the resulting competition leads to an exotic quantum phase tran- 
sition described by a U(l) x Zjv conformal field theory with a fractional central charge. The emerging Zjv 
parafermionic critical degrees of freedom are highly nontrivial in terms of the original atoms or polar molecules 
of the model. We further determine the main physical properties of the quantum critical point in a double-tube 
system which has central charge c = 3/2. In triple-tube systems, we show that the competition between the two 
antagonistic processes is related to the physics of the two-dimensional Z3 chiral Potts model. This work opens 
the possibility to study the exotic properties of the Zjv parafermions in the context of ultracold quantum Bose 
gases. 
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I. INTRODUCTION 

Ultracold atomic systems have become one of the most 
experimentally flexible systems to study strongly correlated 
physics. The possibility of tuning the contact inter-particle 
interactions by varying the scattering length through Fesh- 
bach resonances has allowed the observation of many interest- 
ing phenomena like the BEC-BCS crossover in atomic Fermi 
gas .12 

On top of this wide tunability of the interaction strength, 
the possibility of controlling the shape of interactions through 
the long-range and anisotropic nature of dipole-dipole interac- 
tions has recently attracted considerable attention. In particu- 
lar, dipolar interactions with long-range anisotropic character 
have been observed in chromium atoms by exploiting its large 
magnetic moments.^ An alternative route is to consider het- 
eronuclear polar molecules with large electric dipole moments 
associated with their rotational excitations.^'* The resulting 
long-range interactions between such polar molecules can be 
tuned using dc and ac electric fields.'^ Dipole-dipole interac- 
tions can then be much stronger than the superexchange in- 
teractions between ultracold atoms. This opens an avenue to 
realize a plethora of interesting quantum phases governed by 
long-range interactions. On the one hand, a variety of exotic 
phases have akeady been predicted to occur in dipolar quan- 
tum gases such as topological p-wave fermionic superfluids^ 
and quantum nematic fluids.**^ On the other hand, as shown 
recently in Ref 10 ultracold polar molecules become promis- 
ing candidates for the simulation of condensed matter phe- 
nomena and could provide a robust toolbox for quantum in- 
formation processing. 

In this paper, we investigate the emergence of an exotic 
quantum criticality in coupled one-dimensional (ID) quan- 
tum bosons tubes with dipole-dipole interactions. The stan- 
dard quantum critical behavior of ID quantum bosons is the 
well-known Luttinger universality class.'^^^ Its low-energy 
properties are governed by a relativistic free massless bo- 
son field <i> with a sound velocity v and a Luttinger param- 



eter K which depend on the density and interactions of the 
underlying ID quantum Bose gas. A hallmark of the Lut- 
tinger phenomenology is the power-law decay of the physi- 
cal quantities with non-universal exponents related to the Lut- 
tinger parameter K. Such universality class can in turn be 
viewed as a conformal field theory (CFT) with U(l) global 
continuous symmetry and central c harge c = 1 associated 
to the gapless free-bosonic mode <i).Q32lA more complicate 
quantum-critical behavior has been found recently with ultra- 
cold bosons with the possible realization of the 2D Ising uni- 
versality class.CSEIl In that case, the criticality emerges at a 
quantum critical point characterized by a Z2 CFT with cen- 
tral charge c — 1/2. The resulting gapless degree of freedom 
is a free massless Majoran a (re al) fermion which is half of 
a Dirac (complex) fermion."^'*^' Here, in the context of cou- 
pled ID dipolar bosons, we analyse the stabilization of exotic 
quantum critical points which are described by the SU(2)7v 
CFT with central charge c = 3N/{N + 2)?^ Equivalently, 
the emergent quantum criticality can be captured by the prod- 
uct of a conventional U(l) CFT with the parafermions 
CFT, with central charge c = 2{N - 1)/{N + 2), which is 
known to govern the critic al properties of the 2D Zjv general- 
ization Ising models. Such nontrivial quantum criticality 
for N > 2 cannot be described by a simple free-field theory 
in terms of ID gapless bosons or fermions as it is the case for 
the Luttinger and Ising universality classes. In this respect, 
the resulting critical degrees of freedom are highly nonlocal 
with respect to the original atoms or molecules. Recently, the 
parafermions CFT has attracted a lot of interest in the 
context of non-Abelian fractional quantum Hall states. The 
so-called Read-Rezayi states are described in terms of the Zat 
CFT and their excitations display non-Abelian statistics which 
stem from the parafermionic behavior of these states.l^Sl Our 
work opens the possible realization of this exotic CFT with 
fractional central charge in the context of ultracold atoms or 
polar molecules. 

The experimental setup that allows to create a system of 
coupled ID dipolar bosons tubes, is similar to the one de- 
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FIG. 1: (color online) Planar array of ID tubes of dipolar boson gas 
positioned along the x axis. The dipoles are all aligned along an 
additional external field. 



scribed in Refs l26lj30l A planar array of ID tubes of dipolar 
boson gas, as shown in Fig. [T] is formed by an anisotropic 
2D optical lattice in the z ~ y plane. The strong confine- 
ment in the z-direction completely suppresses the hopping in 
that direction, while we allow a weak hopping between the 
tubes along the y-direction. Neglecting the inhomogeneity of 
the trapping potential, the system is thus a stack of 2D ar- 
rays, containing identical (straight) tubes of bosons, all ori- 
ented along the a;-axis and separated by a distance d = X/2 
in the j/-direction (A being the wavelength of the laser) . The 
dipole moments of the bosons are aligned by an additional ex- 
ternal field which is electric in the case of polar molecules or 
magnetic for chromium atoms. The dipole-dipole interaction 
strength can be tuned by using a rotating field, which induces 
a precession of the dipoles around the z-axis.^ ^' 

The resulting effective Hamiltonian for this system can be 
written as; 



^ ^intra Winter- 



(1) 



^intra IS the Hamiltonian describing the bosons' motion and 
interactions within the same tube: 



n=l 
N 

+ Y. dXidX2 V^f" (Xi - X2) Pn {Xl) Pn (^2) , (2) 
n=l 

where ^fj^ (x) is the bosonic operator that creates a dipole at 
position x in tube n, and pn{x) — 5'Jj(a;)5'„(a;) is the density 



operator In Eq. ([2|, M is the mass of the bosons and VJ^*"^** 
denotes the intratube dipole-dipole interaction. The coupling 
between the tubes is composed of the hopping along the y 
direction, with amplitude J±, and the intertube dipole-dipole 
interactions described by VJJJ'"'^' 



N-l 



n=l 

N~l „ 

+ ^ / dxidX2 V^f"'' {Xi - X2) Pn {Xl) Pn+1 {X2) .(3) 
n=l 

Here, we only consider coupling between nearest-neighbor 
tubes as in Refl26l Unless stated otherwise (in Sec. IV), 
open-boundary conditions in the transverse direction (y) are 
also assumed in Eq. (|3]l. The two terms in Eq. (j3]l are of very 
different nature and can be tuned independently: the depth of 
the optical lattice potential in the y direction controls the inter- 
tube hopping amplitude J±, while VjJJ'^' is controlled by the 
rotating polarizing field. 

In this paper, we will investigate the competition between 
these two sources of coupling, at zero temperature, and in the 
ID regime where the total length in the transverse direction, 
i.e. Nd, is finite. In this respect, small values of N, i.e. N = 2 
or 3, can be obtained in principle, by selective evaporation of 
the tubes by using electric or magnetic field gradients. Alter- 
natively, double-tube system can also be engineered by con- 
sidering double-well optical lattices. ■'^ 

The rest of the paper is organized as follows. In Sec. II, we 
follow a low-energy approach for the model ([T]) in the gen- 
eral N case, using the phenomenological bosonization of the 
bosons" supplemented by CFT techniques. There, we re- 
veal the existence of a quantum critical point in the SU(2)7v 
universality class, which describes the competition between 
hopping and intertube interactions. The properties of this 
fixed point and its stability with respect to general perturba- 
tions are further analysed in the special cases = 2, 3 in Sec. 
Ill and IV. In particular, we show that the low-energy proper- 
ties of triple-tube systems are related to the physics of the Z3 
chiral Potts model. Finally, Sec. V contains our conclud- 
ing remarks. The paper is supplied with an Appendix which 
provides some important technical details. 



II. QUANTUM PHASE TRANSITION IN THE GENERAL 
iV TUBE CASE 

In this section, we will investigate the nature of the quantum 
phase transition which stems from the competition between 
the hopping term and the intertube density-density interaction 
of Eq. ([3]l in the general N case. 



A. Phenomenological bosonization approach 

Let us first specify the form of the dipole-dipole interactions 
which appear in Eqs. (j2] [3]l. The dipole moments are polar- 
ized along a direction d = (sin cos 0, sin6' sin cos 0) in 
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the spherical coordinates (see Fig. [TJl. The energy due to the 
interaction between two dipoles at relative position r is: 



Vdd (r) 



(4) 



where the coupling constant is Cdd — P^/eo for electric 
dipoles of strength p, and Cdd — /^oM'^ for magnetic dipoles 
of strength /i; eq /^o respectively the vacuum permit- 
tivity and permeability. For two dipoles in the same tube, the 
dipole-dipole interaction is then: 

3 sin^ 6 cos^ (h 



intra 
dd 



{Xi - X2) 



Cdd 1 
47r 



\Xi ~ X2\ 



(5) 



It may be convenient to fix 6 to its magic angle value: sin 9 = 
1/a/3 in order to get repulsive interactions. Within this hy- 
pothesis, as in Ref '26', we approximate the dipole-dipole in- 
teraction between two nearest-neighbor tubes by its leading 
contribution: 

T/intcr / \ ^dd COS^ (f> 

Vdd (^1 - ^2) ^ ^ ~d3~ ~ ^2) • (6) 

As we will show thereafter, repulsive intertube dipole-dipole 
interactions are not crucial for the existence of the quantum 
phase transition we investigate. Thus, without loss of gener- 
ality regarding the nature of the intertube interactions, we will 
choose 9 ^ (p = TT /2, i.e. the dipoles are polarized along the 
y-axis, perpendicular to the tubes axis; we then have a short- 
distance intertube attraction 

Trinter / t, ^ ~ ^^t^ 

and a repulsive intratube interaction: 

Cdd 1 



r(5 (Xl - X2) 



(7) 



V^T^X^-X2) 



An \xi - X2\ 



(8) 



With these specifications, we are now in position to in- 
vestigate the physical properties of model ([T]l by means of 
the phenomenological bosonization (or harmonic-fluid ap- 
proach). It is a powerful method that allows to study the low- 
energy properties of strongly interacting ID quantum bosons 
systems.EHi^xhe density operators p„ (x) ,n = 1, . . . ,N can 
be expressed in terms of N bosonic quantum fields (pn{x) 
through;'"* 



Pn (x) 



1 



PO 



00 

E ' 

m— — 00 



(9) 



where are non-universal constants and po is the average 
density in each tube around which the fluctuations take place. 
We introduce the dual bosons associated to the fields ipn 
such that [(pnixi),dx'dm,{x2)] = iSnrnS{xi - X2). They al- 
low us to describe the phase fluctuations of the boson creation 
operator '^l^: 



1 



Po 



■dx<Pn {x) 



1/2 



exp (I- 



(10) 



prn being non-universal constants. In practice, one often only 
needs the leading parts of the expansions ( |9|10| i: 



Pa + —f=dx'Pn {x) 
V 

+2/90 COS ^\/47r [x) + 2itpqx!^ 

'p^ CXp {i\/TT-dn {x)) ■ (11) 



Pn (x) 



*I ix) 



After investigation of the effects of these leading parts, the 
higher harmonics and less relevant terms can then be studied 
a posteriori. 

In the absence of intertube coupling, the low-energy prop- 
erties of the Hamiltonian Q are described by the Tomonaga- 
Luttinger modeli^ 



N 



n 



intra ^ 



= EM- 



K {dxipn 



(12) 



where the velocity v, and the Luttinger parameter K depend 
on the microscopic details of the original model. The Lut- 
tinger parameter has been determined numerically as function 
of the dipole coupling constant Cdd and the average density 
pa-^'^^* In sharp contrast to bosonic gases with contact inter- 
actions where K > 1, K can be much smaller than 1 for 
dipolar interactions. 

The leading-part of the intertube interaction Q can be ex- 



pressed in terms of the boson fields ipn and {)„. using Eq. ( 1 1 



JdxY^l 



9± 



9d 



cos 



^7, 



+ A 



N-l 
dx ^ 
n=l 



/ ^ dx^Pn+ldxV' 



(13) 



a ^ 1/po being the short distance cut off (the average inter- 
particle distance). For the attractive dipole-dipole interaction 
(|7]l, the coupling constants read as follows : 



91- 
9d 

A 



2poa'^J± 

CddPo^ 

Cdd 



(14) 



27r2d3 ■ 

At this point, it is worth noting that the sign of the two co 



sine operators of Eq. ( 13 1 can be changed by the canonical 
transformation on the bosonic fields: 



^2 



2n+l 



^2n+l- (15) 



SO that the case of attractive and repulsive intertube interac- 
tions, such as (j6]) in particular, can be treated on the same 
footing by using the transformation ( 15 i. 
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Let us discuss the different terms that appear in the in- 
teracting Hamiltonian ( [T3] l. They are of very different na- 
ture. The forward scattering term with the coupHng constant 
A will renormalize the velocity and Luttinger parameter of the 
bosonic modes. Now, the structure of the remaining terms 
is intriguing. They describe the competition between the 
bosonic fields iy9„ and their dual fields i9„. The low-energy ef- 



fective Hamiltonian (13 1 is, in fact, similar to the competition 
between superconductivity and char ge-den sity wave (CDW) 
in spin-gap electronic ladder systems,E^126J jj^g ^^^.^^ terms 
in Eq. ( 13 i being respectively the analogues of the Joseph- 



son and CDW interactions. Their scaling dimensions with re- 



spect to the Luttinger-liquid fixed point ( 12 1 are respectively 
A_L = l/(2i^). Ad = 2K. For 1/4 < K < 1, which is a 
regime that can reached with dipole-dipole interactions, they 
are strongly relevant perturbations. The infrared (IR) proper- 
ties of the system is then governed by the effective Hamilto- 



nian: 



He 



n=l •' 

/ dx ^ cos {^/^^ (i?„+i - i?„)) 
^ COS [yArix ((/?„+! - ipn)^ 



(16) 



where the marginal forward scattering term has been dis- 
carded. 

Our strategy thus consists in, first investigating the phys- 
ical properties of model (16i, and then, analysing the main 



effect of less relevant perturbations omitted in (16i, such as 
the forward scattering and higher harmonics of the bosonized 
representations ( 9|10 1. 

One can anticipate on symmetry grounds that the strongly 



relevant terms of model ( 16 1 will not introduce a spectral gap 



for all degrees of freedom. Indeed, it is straightforward to see 
that the model is invariant under two independent global U(l) 
transformations; 



•^n V« + a 

1?« -> +13, 



(17) 



a, f3 being real numbers. Introducing left and right compo- 
nents of the bosonic fields fnL.R — {fn ± '?n)/2, the trans- 



formation ( 17 1 gives a U(1)l x U(l)i^ global continuous sym- 



metry of model ([16]). We thus expect the model to be gapless 
with one bosonic field protected by the symmetry 



17i 



One way to reveal the emergence of this bosonic mode is to 
single it out by performing a change of basis on the bosons 
fields. In this respect, we thus introduce the bosonic field 
^OR,L, and - 1 other fields (I = I, ■ ■ ■ , N - 1) as 

follows; 



1 



OR{L) 



" ^1(^1 _^ ^) ('^1 + • • ■ + ~ 1'Pi+i)r{l) (18) 



Note that in the context of condensed matter physics, the bo- 
son ^oR,L would be called the charge boson since it is simply 
the sum of all bosons and describes fluctuations of the total 
charge. The inverse transformation can be easily found; we 
do not however need its explicit form but only some of its 
general properties; 



fnR(L) 



R{L) 



(19) 



where $ is the — 1 dimensional vector with components 
(; = 1, . . . , TV - 1) and the N vectors v(") (n = 1, . . . , TV) 
satisfy the following relations; 



N 



n=l 



N 







1 

iV 



(20) 



with p,q = I, . . . , N — 1. In the new basis, the effective 



Hamiltonian ( 16 1 simplifies as follows; 



U 



off 



dx 



1 

K 



{dxQo) 



dx 



V 



2 

N-l 



/d. ^[^cos(V^(v(-^)-v(")).0) 

•' n=l 



9d 



(21) 



where 8o and are respectively the dual fields of $o and 
It is interesting to observe that the low-energy effective 
theory ( [2T| i has been found in total different contexts. On 
the one hand, the deconfining phase transition of the 2h-1- 
dimensional SU(iV) Georgi-Glashow model is controlled by 
model (21 1.^*^ On the other hand, as shown recently, it de- 
scribes the Read-Rezayi sequence of non-Abelian fractional 
quantum Hall states.i^Si We explicitly observe, from Eq. (21 



that the $0 field decouples from the interaction terms and re- 
mains gapless. We now need to analyse the fate of the re- 
maining degrees of freedom. The low-energy properties of 
the model depend on the value of K which controls the com- 



petition between the # and fields of Eq. (21 

When 1/2 < X < 1, the most relevant contribution is the 
dual field term with scaling dimension Aj^ = 1/{2K). 
This field is pinned onto the minima of its sine-Gordon model 
in Eq. (21 1 with coupling constant g± > 0; (0) = 0. In con- 
trast, the bosonic $ field is a strongly fluctuating field being 
dual to 0. The resulting phase is a superfluid phase which is 
described by the following leading asymptotics of equal-time 
correlation functions; 

(*U^)*m(0)) ^ a;-i/(2^^^) 
{Pn{x)p^m - pI-^^- (22) 
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In particular, the 27rpo oscillations of the density p„ are short 
ranged. 

When K < 1/2, the most relevant contribution of Eq. (21 
is now the term with the bosonic field with scaling dimen 
sion — 2K. The bosonic field $ is pinned onto the min- 
ima of its sine-Gordon model with coupling constant gd > 0: 
($) = 0. In this phase, the bosonic field is a strongly 
fluctuating field. The resulting phase is a CDW phase, or 
a density-wave phase, where the intratube density are all in 
phase: 



/ / N 1 . COS (27rpoa;) 

(P« (a:) (0)) - Po + ^ 



(23) 



^ being the correlation length which stems from the spectral 
gap opened by the relevant perturbation of Eq. (21 1 with cou- 
pling constant g^, and A is a non-universal constant. While 
the bosonic creation operator has short-ranged correlations, 
this phase allows the formation between the tubes of a molec- 
ular superfluid instability made of iV bosons. Indeed, us- 
ing the bosonization approach 



1 1 1, we find that the operator 



Mt = ^\^\... is described by A/t ^ e*^^"=i = 
We deduce thus the leading asymptotics of molec- 



ular superfluid correlation function: 



(24) 



This instability dominates the density wave (23 i only if i^T > 
iV/2. The situation is in close parallel to the formation of 
molecular fermionic superfluid instability in ID multicompo- 
nent cold fermionic atoms .I^^El 

For decoupled tubes, it is well known that the superfluid 
correlations are dominant when X > 1/2 while the CDW is 
the leading instability for K < 1/2. There is a smooth cross- 
over at A' = 1/2 between these two regions. But as we can 
see from Eqs. ( [22l[23| ), here the situation is in stark contrast 
to the decoupled tubes case, since K = 1/2 marks the onset 
of a zero-temperature quantum phase transition. At K = 1/2, 
the competition between the bosonic fields and their dual is 
maximal since the two strongly relevant perturbations have 
the same scaling dimension Aj^ = Ad = 1: 

N . ^ 

- / E [f^ {f^ (^"+1 - 
n— 1 

+ ^cos (^\/27r((y3„+i - ^„)j . (25) 



B. Non-Abelian symmetry 



The effective Hamiltonian ( 25 1 at K = 1 /2 displays a hid- 

j 1 351361 jjjjg respect, it is 



den SU(2) symmetry when g\^ 



worth introducing the following family of SU(2) matrices: 



1 

71 



(26) 



V 



with n = 1, . . . , iV. The interaction part, V^g, of model (25 1 
can be expressed in terms of these matrices: 



T4ff = / dx 



1 

2^^2 



- (.9-L + gd) Tr (.g.l+i5n) 



(27) 



being the third Pauli matrix. We conclude that the inter- 
action term (27 i displays a non-Abelian symmetry SU(2)i x 
SU(2)i^ for g\^ = gd since it is then invariant under the sym- 
metry: g„ U gnV, U and V being two independent SU(2) 



matrices. When g±^ ^ gd, due to the last term of Eq. (27 1, this 



non-Abelian symmetry is broken down to U(1)l x \5(1)r in 



full agreement with Eq. ( 17 1. 

It is then natural to expect that the quantum phase transi- 
tion between the superfluid and CDW phase occurs precisely 
at g±^ = gd where the effective Hamiltonian (25 1 enjoys an 
enlarged SU(2)i x SU(2)ij global symmetry as well as a self- 
dual symmetry <^„ -H> i?„: 



N 



n* 



n=l 

» N-l 

- % dxY, [cos(V2^(i?„+i-i9„)) 

n+1 H^n) ) 



(28) 



The existence of this global SU(2)i x SU(2)^ continuous 
symmetry leads us to expect that model ( [28] l displays quan- 
tum critical properties which are governed by some CFT Due 
to the SU(2)i X SU(2)/; invariance, the natural candidate is 
the SU(2)fc CFT with central charge c = 3fc/(fc + 2), where 
fc is some integer to be determined. In the next subsection, 
we are going to show non-perturbatively that the model ( 28 1 
is indeed conformally invariant for all sign of gd, and belongs 
to the SU(2)jv universality class. 



C. Conformal embedding approach 

Let us now consider the CFTs underlying the physics of 
our system at the Luttinger-liquid fixed point, i.e. when the 
intertube interactions are turned off, g± = gd = 0. 

At A' = 1/2, each boson field (fi„ {n = 1, . . . , N) is a com- 
pactified boson with radius R = 1/ \/2tt: ^ ipn + ^ttR, 
^ 1^71 + 2ttR . This equivalence can be viewed as gauge 
redundancy in the description since the density and creation 
operators ( 9|10 1 for AT = 1 /2 are invariant under that transfor- 
mation. At this special radius, each boson field describes 
an SU(2)i CFT with central charge c — 1. These CFTs are 
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each generated by left and right SU(2)i currents jnB.,L which 
read as follows in terms of the left-right moving bosons (see 
Appendix Sec. 1 for more details): 



JnR,L 



JnR,L 



1 

27ra 
1 

V2n 



exp 



nR,L- 



(29) 



The SU(2) matrix p6] ) is the spin-1/2 SU(2)i primary operator 
with scaling dimension 1/2. The conformal symmetry of the 
non-interacting part of model (28 i is then S\5(2)i x SU(2)j x 
... X SU(2)]^ with central charge N (i.e., N gapless bosonic 
modes). The N SU(2)i currents ( |29l ) can then be combined to 
form an SU(2)Ar current: 



N 



'■R,L 



,JnR,L- 



(30) 



n=l 



The next step of the approach is to go from these N SU(2)i 
CFTs totheSU(2)jv CFT with central charge 3A^/(iV + 2). In 
this respect, one can use the following conformal embedding 



to investigate the critical properties of model ( 28 1: 



SU(2)i X SU(2)i X • • • X SU(2)i ^ SU(2)jv x Gn, (31) 

where C^at is a discrete CFT with central charge cg„ = N — 
3N/{N+2) = N{N-l)/{N + 2). It is important to note that 
the latter central charge coincides with the sum of the central 
charges of the — 1 first minimal models!^ 



N{N -I) 
{N + 2) 



N+l 

E 

m=2 



1 - 



6 



m {ra + 1) 



(32) 



The tjjv CFT is thus related to the product A^3 x x . . . x 
M.N+1 where A^p denotes the minimal model series with 
central charge Cp = 1 — 6/p(p+l); forp = 3, 4, 5 they respec- 
tively correspond to the Ising, TIM (Tricritical Ising Model), 
and Z3 Potts CFTs.IiSlThe precise identification requires a pro- 
jection?: Qn ^ P{M^ X Ma ... X Mn+i), which has been 
described in Ref 45 

The central point of our analysis is to show that the self-dual 



perturbation of model ( 28 1 



N-l 

Vsd = ^ cos(\/27r(v5„+i-(p„))+cos(y27r(i?n+i-'i9„)), 

n=l 

(33) 

does not depend on the SU(2)Ar CFT. More precisely, we 
show, in Sec. 2 of the Appendix, that the operator Vsd is in- 
deed a singlet under the SU(2)7v CFT and a primary operator 
of the Qm CFT with a scaling dimension equal to one. Hence, 
since it is a strongly relevant perturbation, Vsd introduces a 
spectral gap for the discrete Qjsi degrees of freedom, but leaves 
the SU(2)jv ones intact. We thus conclude that model ( |28] l dis- 
plays critical properties in the SU(2)Ar universality class for 
all signs of gd- 

However, we are not guaranteed that the critical properties 
of the initial model ([T]l are controlled by this SU(2)Ar fixed 



point. One has to investigate the stability of the latter fixed 
point with respect to the marginal forward scattering opera- 
tor of Eq. ( 13 1 and the contributions of the decoupled tubes 



limit such as the higher harmonics of the bosonized descrip- 
tion ( 9|10 l that we have neglected. The marginal forward- 
scattering process breaks the non-abelian SU(2)i x SU(2)/j 
symmetry for gx^ — g^ down to U(1)l x . In this re- 

spect, the symmetry of the IR fi xed point is then U(l) x Zat 
since SU(2)Ar/U(l) - Zn^^ The latter CFT with central 
charge = 2{N— 1) / (A''+2) is called the Z n parafermionic 
CFT and describes the multicritical properties of 2D Z^r 
generalization of Ising models.-^ The self-duality symmetry. 



ipn *H> of model (28 1 corresponds to the Kramers-Wannier 
(KW) self-duality symmetry of Z^r models. In fact, the Zat 
symmetry can be defined directly in terms of the bosons at 
K = 1/2: 



'Pn, 



N 



(34) 



with m = 0, 1, . . . , A^ — 1. Using the representation ( [T0| , 
this Zjv symmetry has a simple interpretation in terms of the 
original bosonic creation operators: 



(35) 



In general, the SU(2)jv and Z^r fixed points are fragile since 
they can be destabilized by several relevant primary operators. 



Here, the presence of the Z n symmetry ( 34 1 does protect the 
system from several relevant perturbations. However, primary 
operators such as the thermal operators e^,, with scaling di- 
mensions Ak = 2k{k + 1)/{N + 2), with fc = 1, . . . [A^/2] 
([Af/2] is the integer part of N/2), are invariant under the Zn 
symmetry and might destabilize the IR fixed point. Fortu- 
nately, some of these operators, such as ei, can be killed by a 
special fine-tuning of the coupling constants by imposing the 
self-duality symmetry. Indeed, under the KW symmetry, one 
has Ek — Cfc, and ei disappears along the self-dual manifold 
the model. As a result, for small values of A^, the SU(2)7v 
or U(l) X Zjv fixed point is likely to be stable thanks to the 
control of the dipole strength and the intertube hopping. In 
contrast, for A^ > 5, a strongly relevant self-dual perturbation 
£2 with scaling dimension 12 /{N + 2) will be generated in 
the Z N sector. The resulting field theory which captures the 
quantum phase transition in this case becomes: 



ScS = 5z„ + X (fx £2 (x) 



(36) 



where 5z„ stands for the action of the Z jv CFT. Model ([36| 
turns out to be an integrable deformation of the Z ^ 
The nature of the phas e tran sition depends on the sign of 
the coupling constant Al ^^ l ^^ l For A < 0, the field theories 
(36 1 are massive and the phase transition is of first-order 
type. For A > it is known that model (36 1 has a massless 



renormalization-group flow onto a Kosterlitz-Thouless (KT) 
U(l) gapless phase with central charge c — 1. In the lat- 
ter case, there is an intermediate gapless phase between the 
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superfluid and CDW phases for N > 5. Adding the contri- 
bution of the gapless mode $o> this phase displays extended 
quantum criticahty with central charge c = 2. Unfortunately, 
our approach does not enable us to fix the sign of the pertur- 
bation of model ([36| and we cannot discriminate between the 
two possibilities for the nature of the transition for > 5. 
In the following, we will investigate the nature of the quan- 
tum phase transition for small values of N, i.e., N ~ 2,3, 
which respectively correspond to double-tube and triple-tube 
systems. 



r-0—o—Q—0- 

I I t t I I I I 

FIG. 2: (color online) The CDW order characterized by Eq. the 
two CDW orders in each tube are in phase. 



III. DOUBLE-TUBE SYSTEM 

In this section, we determine the zero-temperature physi- 
cal properties of the double-tube system which corresponds to 
model ([1) with TV = 2 (see Fig. (2)). 



The harmonic-fluid representation (12 13 i simplifies as fol- 
lows in the N — 2 case: 



n 



N=2 



dx 



dx 



A 



91- 



/^(l?2-^l)) 



+ ^ COS (\/47r (v52 - Vi! 
+ X dx dx^P2dxVi- 



(37) 



The next step is to use the change of basis ( [T8] i: 

1 



72 



(38) 



so that we obtain 



'Hn=2 ^ I dx 



Wo 



Aq 



dx 



— 1 dx 



Ai 



91- 



cos ( v2vr9i 



9d 



cos V87r<i>i 



■(39) 



There is a velocity anisotropy and two different Luttinger pa- 



rameters which stem from the marginal term of Eq. ( 37 1 with 
coupling constant A: 



Ko 



K 



K 

- Xv/K 



(40) 



The effective Hamiltonian ([39| has been first studied in Ref38' 
in the context of a bosonic two-leg ladder for incommensu- 
rate filling. It separates into two commuting pieces: T-Ln=2 ~ 
Hq + Hi (["Ho, "Hi] — 0). The first contribution Hq of Eq. 



( 39 1 has the form of a Tomonaga-Luttinger Hamiltonian; this 
signals that the bosonic field $0 is a gapless degree of free- 
dom. All the non-trivial physics is encoded in Hi which 
describes the competition between the intertube hopping and 
dipole-dipole interactions. As already discussed in Sec. II 
in the general N case, there are two different phases which 
stem from this competition. When 1/2 < Ki < 1, Hi has a 
spectral gap Ai due to the cosine term with the bosonic field 
01. As a result, 6i is pinned onto one of the minima of its 
sine-Gordon model, which we can choose to be (6i) = 
ig± > 0). Thus, at low-energy E ^ Ai, the creation opera- 
tor of the bosons simplifies as follows: "^l^ ^ y/pj^ g^VW^ ©o^ 
while the 27r/9o density operator is short-ranged since it de- 
pends on the strongly fluctuating field $1. The resulting su- 
perfluid phase is then characterized by: 



(p„(x)p,„(0)) ' ^° 



47r2a;2 ' 



(41) 



with n, ni ~ 1, 2. In fact, in this phase, one may also consider 
the symmetric combination of the density with 47rpo oscilla- 
tions in each tube: 



n=l 

^ 2po COS {4:7Tpox + VStt^o) cos (VStt^i) 



(42) 



This operator is naively short-ranged in the phase with (Oi) = 
since it depends on the field $1. However, first-order pertur- 
bation theory with the Hamiltonian Hi cancels the contribu- 
tion of cos (\/87r$i) in p'^^P° . The low-energy description of 
that operator in the superfluid phase thus simplifies as follows: 



(43) 



- 2pa cos i^^TTpdX + \/8^$o) 



We deduce that the two-point correlation function of that op- 
erator has a power-law decay with exponent 4X0- Since 
1/2 < < 1, the decay is much slower than the decays 



of the instabilities of Eq. (41 



The second phase is obtained when Ki < 1/2. The most 
relevant contribution of Eq. ( 39 1 is now the cosine term with 
the bosonic field $1 which is pinned: ($1) = (gd > 0). In 
this phase, the dual field 0i is a strongly fluctuating field. The 
resulting phase is a CDW phase where the two 27rpo CDW on 
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each tube are in phase (see Fig. |2]l: 

{Pnix)p^{0)) - p^ + A^^'g^, (44) 

A being a non-universal constant and ^ ^ ui/ Ai is the corre- 
lation length which stems from the spectral gap of the Hamil- 
tonian?^!. In this phase, the dimer instability Af^^ = 'I'|5'2 ^ 
Pf) e^^^ ^" has a power-law correlation function with a larger 
decay 1/Kq (Kq < 1). 

The quantum phase transition between these two different 
phases occurs at Ki = 1/2. From Eq. (40 1, we observe 



that since A < 0, we have Ki < K, i.e., attractive intertube 
interactions such as in Eq. (j7]l are favorable in order to reach 
the transition. The effective Hamiltonian which governs this 
transition writes: 



'Heft - 

dx 



dx 



5-L 



cos 



9d 



cos 



(45) 



This model is well known (see e.g. Ref 49) and can be 
exactly diagonalized by introducing two Majorana fermions 
i}-^ . This procedure is nothing but the standard bosoniza- 
tion of two Ising models. ^ ' The refermionization rules are 
given by 



1 



Tra 
1 



cxp 47r$i/; 

exp ^— i\/47r$iL 



(46) 



where the anticommutation between the right and left-moving 
Majorana fermions is taken into account by the prescription: 
— i/4. The effective Hamiltonian can then be 
refermionized: 



2 



dx [in>,iCUi + «"^2Ciei] , (47) 



where the two masses are given by: 



mi 

7712 



{9d 



9±) /a 
91-) /a- 



(48) 



The effective model (45 i at if i = 1/2 is thus mapped onto 
two decoupled ID Ising models in a transverse field. The 
Majorana fermion ^ always has a positive mass m2 which 
means that, in our convention, the corresponding Ising model 
belongs to its disordered phase where the disorder operator 
condenses: (/i2) 7^ 0. In contrast, the Majorana fermion 
has a smaller mass mi which can change its sign; by tuning 
the hopping g^, we are thus able to reach the Ising Z2 quan- 
tum critical point where mi = 0. 



The different order parameters defined in Eqs. (41 44 1 can 
be expressed in terms of the Ising order (cri.2) and disorder 
(/Xi 2) parameters of the two underlying Ising models. In this 
respect, we need the Ising description of all the vertex opera- 
tors with a scaling dimension equal to 1 /4p3Ell 

cos(^/7r$l) ^ fJ-ifJ-2, y2cos(^/7^9l) 0-1^2 
V2sin(y7f$i) - 0-1(72, \/2 sin ( ^/ttG 1 ^1(72, (49) 

where the Ising order parameters are normalized at the Z2 
quantum critical point according to: 



{aa (r, x) ab (0,0)) 



Sab 



1/8 ■ 



(50) 



with a,b ~ 1,2, and similarly for the two-point correlations of 
the Ising disorder operators. At low-energy, i.e., E ^ TO2, we 
can average the operators of the second Ising model {{^2) 7^ 
and ((T2) — 0). The bosonic creation and density operators 



then simplify as follows, using Eq. (49 1: 




Acqs 



(51) 



27rpo2^ + y27rKo<I>c 



Ml, 



with Kq — 1/2. In particular, we observe that the bosonic 
creation operator expresses in terms of the Ising order operator 
CTi, in full agreement with the identification of the discrete Zn 



symmetry (35 1 made in Sec. II. 

It is then straightforward to deduce the zero-temperature 
phase diagram of model (45 1. When g± < g^t, i.e. mi > 0, 
the Ising model belongs to its disorder phase with (cri) — 
and (/ii) 7^ . We recover the CDW phase with the leading 
asymptotics ( |44| . When g±^ > gd, the Ising model belongs to 
its ordered phase with (cti) 7^ and the double-tube system 
sits in its superfluid phase (41 1. The quantum phase transition 



between these two phases occurs for g± = gd and is driven 
by the Ising degrees of freedom. The resulting critical prop- 



erties are determined by the low-energy description (51 1 and 
correlations (fSOli: 



{^i (x) (0)) 

{Pn {x) Pm (0)) 



-l/(4Ko)-l/4 



pI 



Att'^x-' 



+ A 



cos (27rpo2;) 

rKo + l/A 



(52) 



with Kq ~ 1/2. The quantum phase transition is thus de- 
scribed by a U(l) X Z2 CFT with a central charge c = 
1 + 1/2 = 3/2. The SU(2)2 critical properties of the per- 
turbation (33 I for = 2 are weakly broken down to U(l) x 
Z2 since the bosonic fields $0 and $1 have different Luttinger 



parameters (40 1 and velocities. 



IV. TRIPLE-TUBE SYSTEMS 

In this section, we determine the zero-temperature physi- 
cal properties of the triple-tube system which corresponds to 
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model ([T]l with = 3 (see Fig. [3]). We consider here two dif- 
ferent geometries: one with open-boundary conditions (Fig. 
[3|a)) and the other with periodic boundary conditions, where 
the tubes form an equilateral triangle (see Fig. |3jb)). In the 
latter case, we assume a non-frustrated situation where all in- 
tertube density-density interactions are attractive. 



these degrees of freedom is given by: 

ax 



cos (Vtt 9i) cos ( vStt 62 



cos (-\/7r<&i) cos ( \J 37r$ 



(55) 



The bosonic field $1 is at the free-fermion point, i.e. the terms 
of the form cos or cos (/? 0i) that appear in the Hamil- 

tonian all have /? = -y/7r. For this reason, we can introduce two 
effective Ising models with order parameters u\,2 and disor- 
der fields /ii_2, so as to refermionize the vertex operators of 
Eq. (55 I with a scaling dimension 1/4, just as we did for the 
double-tube system (see Eq. (|49])): 



(a) 



(b) 



FIG. 3: A triple-tube system for (a) open boundary conditions; (b) 
periodic boundary conditions. As in Fig. [T| the tubes are separated 
by a distance d. 



A. Z3 emerging quantum criticality 

The harmonic-fluid representation for the triple-tube sys- 
tem of Fig. |3ja) is given by model (12 13 1 with = 3. Here, 
we follow the general strategy of Sec. II and first consider the 
most relevant contribution, i.e. Eq. (I61 with = 3. We 



single out the gapless U(l) $0 bosonic field by switching to 
the basis ([TSj with iV = 3. The inverse transformation ([T9]l 
explicitly reads as follows: 







1 


(^1 = 




71 






2 


<P2 = 




7i 






1 


^■i = 




71 



$1 



$2 



$1 



1 

71 



$2 



=$2 



with similar definitions for the dual fields i?i.2 3- The effective 



Hamiltonian ( 16 1 for = 3 separates into two commuting 
pieces: T-L^s = Ho + Hs, where Hq takes the form of the 
Tomonaga-Luttinger model: 



■He 



dx — 
2 



A 



(54) 



Us 



IV 



dx 



V2 



dx 



(9, $2)' + (9,62)' 

.g_L0'i/i2 : cos {y^TT 62 j : 



5dA'iM2 : cos 



(56) 



The next step of the approach is to switch to a different ba- 
sis to determine the main effect of the relevant perturbation of 
model (56 1. In this respect, we exploit the fact that the CFT 



defined by the product of one Ising CFT (the one with disorder 
operator /ii) and the CFT associated to the boson $2, can be 
described in terms of the product A^4 x A^5 CFT, i.e., TIM 
X Z3 Potts CFTs. Indeed, the TIM and the Z3 Potts CFTs 
have central charges c = 7/10 and c = 4/5 respectively, so 
that the sum gives c = 3/2, the total central charge of an Ising 
CFT c = 1/2 and a CFT associated with one boson c = 1. 
The TIM CFT describes the physical properties of the two- 
dimensional dilute Ising model at its tricritical point. ''^ The 
precise identification of the conformal embedding has been 
derived by the authors of Ref,45 and requires a projection 
which restricts the primary operators of A^4 x M.^ to the sub- 



(53) set: where H < r < p - 1,1 < s < p) 

denotes the primary fields of the minimal model series M.p 

.Elln 



(P) _ [(p+l)r- 



CFT with scaling dimensions Ar,s — 2p{p+i) 

Sec. 3 of the Appendix, we show that we have the following 
correspondences: 



^1 : cos {yini^-^ : +cri : cos (^3^02^ : CTtj^ (57) 
^1 : cos (|\/37r$2^ : -^i : cos (^V37r82^ : ^ CTTiMeZg, 



The remaining degrees of freedom in "Hs are strongly coupled 
and control the quantum phase transition between the super- 
fluid phase (|22| and the CDW phase ^ at K ^ 1 /2. When 
_ftr = 1/2 and in the new basis, the effective Hamiltonian for mal operator of the Z3 Potts CFT with scaling dimension 4/5 



cttim and crrj^jj^ being respectively the magnetization and sub- 
leading magnetization operators of the TIM CFT with scaling 
dimensions 3/40 and 7/8.'^ In Eq. (57 1, e^., refers to the ther- 
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With all these results, one can now express the low-energy 
effective Hamiltonian in terms of these new degrees of free- 
dom: 



2a 



dx 



{gi- + 9d) M2C^TIM 



+ (~5-L + gd) M20-TIMeZ3 



where -U^^^ and 



(58) 



respectively denote the Hamiltonians 
of the TIM and Z3 CFTs. The operator /i20'TiM ^ strongly 
relevant perturbation with scaling dimension 1 which couples 
the Ising model associated to the Majorana fermion with 
the TIM CFT. The Ising critical point is destabilized upon 
switching on the disorder field ^2, while the dr^jj^ perturba- 
tion is known to be a massive integrable deformation of the 
TIM CFT for all signs of its coupling constantPl We deduce 
that a spectral gap A for the Ising and TIM degrees of freedom 
is opened as soon as the tubes are coupled. In the low-energy 
limit where £^ <C A, the low-energy effective Hamiltonian 



(58 I then simplifies as 



Us = -Ho' +A dx (-.g^ + gd) 



(59) 



A being a non-universal constant. This model is a massive 
and integrable perturbation of the Z3 CFT.I^ A spectral gap is 
thus generated, except at the fine-tuned point = gd where 
the model displays a quantum critical behavior in the Z3 Potts 
universality class. The position of the critical point, i.e. g± = 



gd, corresponds to the self-duality symmetry of model (55 1: 

$1,2 ^ Oi,2- 

Adding the contribution of the gapless boson $o> we de- 
duce that the quantum phase transition of the triple-tube sys- 
tem is governed by the SU(2)3 CFT with central charge c — 
1 + 4/5 = 9/5; this is in full agreement with the general 
analysis of Sec. II. 

As in the N = 2 case, the marginal forward- scattering op- 
erator Ofs — dxip2{dxVi + dxfTi) is expected to introduce 
anisotropics and the symmetry of the quantum critical point 
is U(l) X Z3. In particular, using the new basis (53 1, we find 
that 



(60) 



As expected, the first term of this equation merely intro- 
duces velocity anisotropics and different Luttinger parame- 
ters for the bosons $0,1,2 as in the double-tube case. On 
the other hand, the last term gives a residual coupling be- 
tween the gapless boson $0 and the remaining degrees of free- 
dom. Such term has also be encountered in the analysis of the 
zero-temperature phase diagram of ID three-component cold 
fermions; its main effect is not clear.**^ However, we think that 
the latter perturbation will not destroy the quantum criticality 
of the phase transition of the triple-tube model. 

In any case, one can eliminate this term by considering a 
triple-tube system with transverse periodic boundary condi- 



tions as depicted in Fig. |3jb). In that case, the forward- 
scattering perturbation of this problem is: 



Cfi = dxf2{dxfi + dx^ps) + dx(pidxip3 
= {d.^of - \ [{dx^i? + {dx<^2f_ 



(61) 



so that this contribution is exhausted by velocities and Lut- 
tinger parameters renormalization. At the self-dual point 
.9J- = gd, the effective Hamiltonian which controls the transi- 
tion of the model has a similar structure as in the open bound- 



ary case (see Eq. (58 1): 



Hi 



IV 

V2gd 



I dx {endxa - iid^eL) + nr"" + 

2Trgd 



dx ULoCF' 



M20'tIM 



dxieneL- (62) 



We observe that here, the spectral gap for the Ising and TIM 
degrees of freedom is now explicit, due to the mass term 
(m2 — 2TTgd/a) for the Majorana fermion Indeed, since 
this mass is positive, the corresponding Ising model belongs 
to its disordered phase with {112) 7^ 0. Averaging out these 
degrees of freedom, one is left with a TIM CFT perturbed 
by the subleading magnetization operator crrpjj^, which gives 
a spectral gap.l^ We conclude, as in the open case, that at 
low-energy, the Z3 degrees of freedom remain untouched and 
display quantum criticality with central charge c = 4/5. 

Our approach allows us to determine the leading asymp- 
totics of the correlation functions at the quantum critical point 
with U(l) X Z3 symmetry. The bosonic creation operators ^fj^ 
can be expressed in terms of the bosons (53 1 with K = 1/2: 



1,3 



(63) 



with A'o — 1/2. Using Eq. (49i and taking into account that 
the Ising model with disorder parameter /i2 sits in its disor- 
dered phase, we have 



t 

1,3 



t ^ pi^/-^/3Ko Bo iy/4TT/ 3 62 



(64) 



The next step is to express the operators trie* v '^^'^ ®^ and 
g*A/47r/3e2^ with respectively scaling dimensions 5/24 and 
1 /3, in terms of the TIM x Z3 degrees of freedom. We find 
the following identification: 



(65) 



aieV-/3e2 _ aTiMCTZ3 



,i^47r/3 02 



fiTIMO'Za 



CTZ3 being the Z3 spin operator with scaling dimension 2/15, 
and etim, the thermal operator of the TIM CFT with scal- 
ing dimension 1/5. Finally, at low -energy, i.e. for energies 
smaller than the spectral gap of the TIM degrees of freedom. 
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the bosonic creation operators of the triple-tube system 
simplify as follows: 



(66) 



We observe that '^l,^ expresses directly in terms the Z3 spin 
operator in full agreement with the identification ( [35] l in the 
general case. Similarly, we can find the low-energy represen- 
tation of the 2TTpa CDW on each tube pl'^P" at K = 1/2: 



^27rpo 



J2Trpax+iy/4TrKo/3 *o 



(67) 



with Kq ~ 1/2 and /xz., is the Z3 disorder spin operator which 
is dual to the Z3 spin operator . Using the results ( 66 67 1, 



we deduce the leading asymptotics of the equal-time correla- 
tions at the quantum critical point of the triple-tube model: 



(*U^)*m(0)) - X 
{Pn {x) p,n (0)) ^ p\ 



l/(6A'o)-4/15 



cos (27rpoa;) 



with Ko ~ 1/2. The quantum phase transition is thus de- 
scribed by a U( 1 ) X Z3 CFT with central charge c = 1+4/5 = 
9/5. 



B. Stability of the quantum critical point and the Z3 chiral 
clock model 



The next step of the approach is to investigate the stability 
of this quantum critical point to generic perturbations allowed 
by the symmetries of the model. This includes all irrelevant 
perturbations, that we have neglected in the continuum de- 
scription, which may change their status at the new U(l) x 
Z3 fixed point. The entire content of the Z3 Potts model is 
known'^ and the thermal operator e^., is the only relevant op- 
erator having zero conformal spin and preserving the global 
Z3 symmetry of the model. The presence of two independent 
coupling constants in the model, namely J± and Cdd, should 
be enough to kill this relevant operator and reach the Z3 criti- 
cal point. The other relevant operators of the Z3 Potts model 
carry a non-zero conformal spin: $(2/5.7/5) and $(7/5 2/5) are 
Z3 operators with scaling dimension 9/5 and conformal spin 
S = ±1. These non-Lorentz invariant perturbations, if gen- 
erated, might drive the system to the fixed point of the chiral 
three-state Potts universality clasJ^, the general low-energy 
effective field theory of triple-tube systems, which controls 
the vicinity of the Z3 fixed point, would then read as follows: 



-Hi 



nil 

tin 



f^a;5$(2/5,7/5) +'5*$ 



(7/5,2/5) 



(69) 

where 6 and Sg ^ A {—g± + gd) are small coupling constants 
describing the departure from the Z3 quantum critical point. 
It is thus important to find out if the relevant non-scalar op- 
erators $(2/5.7/5), *i'(7/5, 2/5) are permitted in the triple-tube 
systems investigated here. 



In this respect, we observe that a non-zero conformal spin 
perturbation emerges in the continuum limit of the triple-tube 
system with open boundary conditions (Fig. 3(a)), when ex- 
panding the square root in the harmonic -fluid representation 
of the bosonic operator ( 10 1; for K — 1/2, it reads: 



Os=±i = ^ (d^^n+i + d^iPn) COS (^\/27r (??„+i - 

(70) 

This term is a marginal perturbation with conformal spin S — 
±1, and it is invariant under parity, since x — >■ —x , ipn{x) — > 
—tp„ {—x), which is indeed a symmetry of Eqs. (9 10 1. In 



the triple-tube model with periodic boundary conditions (Fig. 
3(b)), the non-Lorentz invariant perturbation becomes: 



s=± 



1 = ^ {d^ip„+i + do^ipn) COS (^27: (??„+i - i9„)^ 



(71) 

Using the basis (53 1 and the refermionization (46 49 1, we fo- 
cus on the following combinations for each triple-tube system: 



1 = 4=M2Cri : 5a:$2 COS (^37102^ 



(1) 

s=±i ^ 



(72) 



where the second operator occurs for the triple-tube model 
with periodic boundary conditions. We can now express these 
non-Lorentz invariant perturbations in terms of the TIM x Z3 
degrees of freedom. We find the following identification up to 
some unimportant normalization factors: 



O 



(1) 

s=±i 



Ai2CTTIM [$(2/5,7/5) + $(7/5,2/5)] 



9x$2iCi?,Ci ^ CTIM [$(2 



/5,7/5) 



$ 



(7/5,2/5)J 



(73) 



Thus, at energies much smaller the spectral gap of the Ising 
and TIM degrees of freedom, the physical properties of triple- 



tube systems are governed by the effective Hamiltonian ( 69 1 
with 5^6* (and 6g = A{—g± + gd))- The latter Hamiltonian 
is known to describe the physical properties of the 2D uniaxial 
Z3 chiral Potts model, the so-called Ostlund-Huse model,^^ 
in the vicinity of the three-state Potts critical point. Inter- 
estingly enough, the competition between superfluidity and 
CDW in triple-tube systems is thus sim ilar to the domain-wall 
wetting transition in chiral clock model.l^We note that a ID 
hardcore boson chain also has physical properties described 
by the Z3 chiral Potts model.Esl 

The Ostlund-Huse model on the square lattice is defined as 
follows: 



n 



OH 



cos 



27r 



{rii — rij + A) 



<i,k> 



27r 



{Ui - Uk) 



(74) 



where a variable rii = 0, 1, 2 is associated with each site i of 
a square lattice and the sum < i,j > (respectively < i,k >) 
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1/2 



FIG. 4: (color online) Schematic phase diagram of the 2D Z3 chiral 
clock model; C, L, and I denote respectively the commensurate, liq- 
uid, and incommensurate phase; Tc is the critical temperature of the 
three-state Potts model and A the asymmetry parameter. 



are taken over nearest neighbours in the axial, e.g. x, (re- 
spectively y) direction. The Ostlund-Huse model displays a 
global Z3 symmetry and has been introduced to describe the 
commensurate-incommensurate transition observed in mono- 
layers absorbed on rectangular substrates. For A = 0, the 
model reduces to the ordinary 2D Z3 clock model, which 
is equivalent to the three-state Potts model, and enjoys a S3 
(the permutation group of three objects) global symmetry. 
Such a model displays a second-order phase transition from 
a commensurate (or ferromagnetic) phase to a liquid (or para- 
magnetic) phase in the Z3 parafermionic universality class. 
The phase diagram of the model with uniaxial anisotropy 
A 7^ has proved complicated and very controversial (see 
Refs'59'60' for a review). The main effect of the asymme- 
try parameter A is to introduce a new incommensurate gap- 
less phase (the so-called floating phase, with central charge 
c — 1) between the conventional commensurate and liquid 
phases (see Fig. 4 for a schematic phase diagram). The 
floating phase is characterized by a modulated order with a 
wavevector which varies continuously with T and A. There 
is considerable controversy over the possible existence of the 
Lifshitz point (point B in Fig. 4). To the best of our knowl- 
edge, the situation is still not settled. Three different scenari 
are possible: 

(a) There is a Lifshitz point at finite value of A and the Potts 
critica l point (point A in Fig. 4) extends along the AB line of 

pjg ^ 1331611621 

(b) The Lifshitz point moves to A = 0: there is no direct 
transition between the commensurate and liquid phase but an 
intermediate floating phase.'^^Mli There are two phase transi- 
tions to reach the liquid phase from the commensurate one: 
first, a commensurate-incommensurate transition followed by 
a KT transition. 

(c) There is a Lifshitz point at finite value of A but the AB 



line of Fig. 4 is a new critical line with critical e xpon ents 
different from the three-state Potts universality class.ESSl 

In our triple-tube systems, the case (a) corresponds to a 
U(l) XZ3 quantum phase transition between the CDW and 
superfluid phases, described by the properties ( [68] l. In the 
second case (b), there is no direct transition between these 
two phases but an intermediate incommensurate gapless phase 
will then emerge. Adding the contribution of the gapless bo- 
son this incommensurate gapless phase has central charge 
c = 1 + 1 = 2. The correlation functions ( |68| l, in this float- 
ing phase, will be non-universal and incommensurate with a 
wave-vector Q ^ 2-kpq. We expect the extent of this incom- 
mensurate phase to be tiny since 5 = d* is, small in our ap- 
proach. Finally, in the last case (c), we still have a direct quan- 
tum phase transition between the CDW and superfluid phases. 
The correlations functions, at the quantum critical point, are 



still described by Eqs. (68 1 but with different (unknown) uni- 
versal exponents. 

Unfortunately, the low-energy properties of the effective 
field theory (69 1 with 5 — 5* and 5a = Q are not known and 



cannot help to resolve the controversy.*^ However, on gen- 
eral grounds, a relevant 5 = ±1 conformal spin perturbation 
usually gives rise to some incommensurate behavior.i^^^El 
The second case, i.e., the existence of an intermediate incom- 
mensurate gapless phase, is thus maybe the most probable. 
In addition, we expect a different behavior between the two 
triple-tube systems of Fig. 3 since they do not share the same 
symmetry. In the equilateral triangle situation, we have a lat- 
tice global S3 invariance in contrast to the open case (Fig. 
3(a)). Since the Ostlund-Huse model (|74]i does not have this 
invariance except at A = 0, it is natural to conjecture that 



5 = 5* = in the low-energy effective Hamiltonian ( 69 1 for 
the periodic case. In that case, the quantum phase transition 
between the superfluid and CDW phases belongs to the U(l) 
XZ3 universality class. Unfortunately, within our approach, 
we cannot directly check that 5 = 5* = Q. Indeed, the evalua- 
tion of these coupling constants requires the exact knowledge 
of the non-universal amplitudes on Ising and TIM degrees of 
freedom in Eq. (73i for the theory ([62]l. In this respect, it 



would be very interesting in the future, to numerically investi- 
gate the zero-temperature phase diagram of lattice triple-tube 
models of Fig. 3. In particular, it would further shed light on 
the nature of the quantum phase transition and its relation to 
the physics of the 2D Z3 chiral Potts model. 



V. CONCLUDING REMARKS 

In summary, we have investigated the low-energy proper- 
ties of a system of N ID dipolar bosons tubes coupled to- 
gether via intertube hopping and dipole-dipole interactions. 
Here, we have focused our analysis on the quantum-phase 
transition for incommensurate filling that reflects the compe- 
tition between the two different coupling mechanisms using a 
phenomenological bosonization approach. The transition sep- 
arates a superfluid phase from a CDW phase when the Lut- 
tinger parameter of the bosons fields is K = 1/2. This value 
can be reached here thanks to the long-range nature of the 
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dipole-dipole interaction. In stark contrast to a single-tube 
case where the transition is a simple cross-over, the main ef- 
fect of coupling the ID tubes through hopping and density 
interactions is to give rise to a genuine quantum-phase transi- 
tion. 

Using the powerful machinery of the CFT approach, we 
have determined the main long-wavelength properties of this 
transition which turns out to be very exotic. In particular, we 
have revealed that the quantum phase transition is described 
by the SU(2)Ar CFT, or equivalently to the U(l) x CFT, 
with a fractional central charge c = 3N/{N + 2). The 
critical degrees of freedom are highly nontrivial and nonlocal 
with respect to the original atoms or polar molecules of the 
model. The quantum critical point can be attained by a fine- 
tuning of the coupling constants that seem reasonably realistic 
in the context of cold polar molecules. In this respect, this 
work opens the possibility to investigate the exotic physics of 
parafemiions in the context of ultracold quantum bosonic 
gases. 

In the N = 2 case, we have investigated in details the 
main characteristics of the quantum critical point which was 
shown to belong to the SU(2)2 universality class with cen- 
tral charge c = 3/2. The triple-tube case (N = 3) is par- 
ticularly promising since the quantum phase transition is gov- 
erned by U(l)xZ3 degrees of freedom. In particular, within 
our low-energy approach, we have connected this problem to 
the physics of the 2D Z3 chiral Potts model. In this respect, 
it will be important to carry out a thorough large-scale nu- 
merical analysis for mapping out the zero-temperature phase 
diagram of the triple-tube systems by adding an optical lattice 
in the x-direction. The commensurate case is also interesting 
and will be investigated elsewhere. We hope that ongoing ex- 
perimental research on polar molecules will allow to probe the 
quantum phase transition discussed in this paper. 



c = 1, generated by the left and right SU(2)i current: 



JnR,L 



JnRX 



1 

1 



: exp(T«v87r ipnR,L) ■ 

dx'PnR,L, 



(Al) 



where : ^ : is the standard normal ordering of a bosonic op- 
erator A. In this work, the chiral bosonic fields ^PnR.(L) ^6 
defined as: 



1?. 



^nL + VnR, 
VnL — 'firiR, 



(A2) 



and we are working with the prescription: [ipnR, fmL] = 
'iSnm/4:. The bosonic fields are normalized by the following 
operator product expansion (OPE): 



VnR{z)(PmRil^) 



in 



ln(z — uj), 



(A3) 



with z — VT + ix, z — VT — ix, and t is the imaginary time. 

Eqs. (jATJ give a free-boson representation of the SU(2)i 
currents that satisfy the so-called SU(2)i Kac-Moody algebra 
defined by the OPeP 



87r2(z-a;)2 ^ 27r(z - ' 



(A4) 



with a, /3, 7 = X, y, z and a similar result for the right cur- 
rent. Indeed, using Eq. ( |A3| l and OPEs involving vertex 



operators 1 



: e 



{z):e 



ibipnL 
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Appendix A: Free-field representations 

In this Appendix, we present some important technical de- 
tails on the SU(2) n quantum critical properties of the effective 



Hamiltonian ( 25 1, which controls the quantum phase transi- 



tion of the coupled tubes system. 



1. Free-field representation of SU(2) jv CFT 

As stated in Sec. II B, for g± = gd and K = 1/2, the 
Hamiltonian is given by Eq. (28 1 and displays an enlarged 
SU(2)i X SU(2)fl global symmetry. Let us consider its non- 
interacting part where denotes the g± = point. 
Each boson ipn describes an SU(2)i CFT with central charge 



one shows that: 

JnL(^)j«L(w) 



Att{z — Ul) 



(A5) 



1 



4tt^{z - w)2 

1 



± 



1 



47r2(z — w)2 tt{z — uj) 

1 tV&tt 
27r^/27^47^(z-w) 

1 

87r2(^-a;)2- 



(A6) 

This is nothing else but the algebra ( |A4[ ). 

As any CFT, the SU(2)i CFT described by the bosonic 
fields fnR(L) possesses a stress-energy tensor. Its left - i.e. 
holomorphic - component is given by: 



47r2 



■ JnL ' JnL ■ ; 



(A7) 
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and satisfies the defining relation of a stress-energy tensor of 
aCFT:'^ 

rWTM^^ + ^ + ^^, (A8, 

where d — d/dz, and c is the central charge for the CFT in 
question; for the SU(2)i CFT, c = 1 as already mentioned be- 



fore. Using the representation ( Al i, the SU(2)i stress-energy 
tensor can be directly expressed in terms of the boson field: 



r„ = 27r : {d^ipubY 



(A9) 



it is straightforward to check that the definition ( |A8| l is indeed 
reproduced. 



The non-interacting part of the Hamiltonian (28 1, "Hq, can 
be written in terms of the SU(2)i currents ( |A1[ ): 



N 



/2Trv 
dx — (: j„L • j„L : + : • Ur ■) , (AlO) 



and its underlying CFT is thus SU(2)ixSU(2)iX---xSU(2)i 
with central charge N . If we now combine the N SU(2)i 
currents together: 



1 ^ 

— : exp(=Fi\/87rV'nfl,L) : 

ZTT ^ — ' 



n=l 
N 



(All) 



the resulting object is an SU(2)Ar current. Indeed, one can 
show, similarly to the SU(2)i case, that it satisfies the SU(2)Ar 
Kac-Moody algebra: 



87r2(z-Lj)^ 27r(z-w) ' 



(A12) 



The holomorphic part of the stress-energy tensor of the 
SU(2)^ CFT is defined by: 



A^ + 2 



It I 



L ■ 



(A13) 



It satisfies the relation (AS i with a central charge csu(2)jv = 
^N/{N + 2). Using theldentification (|A11||, we find a free- 



field representation of the SU(2)jv stress-energy tensor in 
terms of the bosonic fields (pnL- 



67r 



N 



'SU(2)j, 



iV + 2 
1 

iV + 2 

27r 

iV + 2 



n=l 

N 



: COS [\/87r((^„L - (/s^l)] : 



n^m— 1 
N 



Y ■ d^'^nhdx'-PraL ■ ■ (A14) 



2. Self -dual perturbation and Qn primary field 



In Sec. lie in order to investigate the critical properties 



of the model, we have introduced the Qn CFT through the 



conformal embedding ( 3 1 



SU(2)i X SU(2)i X • • • X SU(2)i ^ SU(2)jv x Gn- (A15) 

Its stress-energy tensor is the difference between the 
SU(2)i and the SU(2)7v stress-energy tensors, i.e., in the 
bosonic description: 



N 



SU(2) 



(A16) 



2tt{N- 1) 



iV + 2 
1 



N + 2 
N + 2 



: COS (V87r(</?„L - <PniL)) ■ 



: dxlfinLdxiPniL 



(A17) 



n^rn 



Let us show that indeed Tg„ satisfies the OPE of stress-energy 
tensors ([A8|) with a central charge cg„ = N{N ~1)/{N + 2). 



Since we already know the OPE that T„ and Tgjj(^2)N satisfy 
(see Eq. ( |A8[ l), we only need to compute the OPE between 
the N SU(2)i and the SU(2)Ar stress-energy tensors. Using 
the free-field representations ( |A14 A17i and Eq. (A5i, one 
obtains: 



N 



Tsv(2)„{z)Yt„{uj) 



3N 



1 



2T. 



SU(2)^ 



The OPE ^n(^)^su(2)jv (i^) can be easily deduced from 

the former by exchanging z and a;, then performing an expan- 
sion of z around w; it has the exact same form. Combining all 
these OPEs together, we finally have: 



2(7V + 2) {z-uY 



(z -w)2 
(A18) 



Tq^{z)Tq^{uj) 



N 



n,m— 1 



N 



N 



SU{2) 



-7SU(2)„(2) Y ^"(^) ~ E ^»(^)^SU(2)„(w) 



N{N~1) , 2Tg,iLu) , dTg,{cu) 



2(7V + 2)z4 (z-w)2 



UJ 



(A19) 



which proves the value of the central charge cg„ = N{N — 
l)/(iV + 2) asserted above. 

In the analysis of Sec. II C, we have used some crucial 



properties of the self-dual term (33 i that is responsible for the 
quantum phase transition in the dipolar bosons tubes. With the 
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definitions above at hand, we now demonstrate these proper- 
ties. In terms of the bosonic fields, the self-dual term reads: 



N-l 



Vsd = ^ : cos(\/27r((^i+i - ipi)) 



i=l 



+ : cos(V27r(t?,+i - 



(A20) 



Let us check explicitely that it is a singlet under the SU(2)Ar 
CFT: 



0, 



(A21) 



and similarly for the right component. First of all, we have, 
using Eq. (|A5[): 



(A22) 



^ 7 ^(^'"''^^[C'k+iCk - Sk+iSk] 

At: ^-^ (z — uj) \ 

i ^ ' 

+<5''''^[Cfc+iCfc-^fe+i5fc])(^,w), 
^ -. — r ( Cfc+iCfc — Sk+iSk 

4:TT[Z — U!) \ 

+Ck+iCk — Sk+iSk^ (i^,'^), 

: cos(y2^(i?fe+i - : 

^ — ^ (Ck+iCk — Sk+iSk 

47r(z — Lu) \ 

+C'k-i-iC'k — Sk+iSkj (i^, w), 
where we have adopted the shorthand notation: 

: cos (Viiripk) Ck, : cos {V27T-dk) := Ck, (A23) 

and with similar definitions for the sine functions Sk and 5^. 
We conclude thus: II^ {z)Vsd{<^,^) ^ 0, and similarly for the 
y component of the SU(2) jv current For the z component, 
we have, using Eq. ( |A5| i: 



iL^'iz) : cos V27r(.^fe+i - ipk) : (w,a;) 



sin %/27r(i^fe+i ~ ipk) : ) (w, 



(A24) 

The same also applies for the second term of Vsd and that 



finally proves Eq. ( A21 



The second important property is that V^d is a primary field 
under the CFT, with scaling dimension Ay = 1: 

TgAz)V.d{u^M - + . (A25) 



Let us check that Vsd displays such property. First, we have 
just shown that lL{z)Vsd{^^,^) ^ 0, so it is straightforward 
to see that T'su(2)„ {z)Vsd (w, w) 0. We are left with the OPE 
of the N SU(2)i stress-energy tensors with Vsd: 



N 



^r„(z) : cos[v/27r(.^fe - Lpk+i)] : {oj,uj) 



N 

E 



1 1 



{^n,k + ^n,k + l) 



n=l ^ ' 

X : [CkCk+i + SkSk+i] ■■ 
1 f 

+ V <^ <5„,fc [diCk)Ck+i + diSk)Sk+i 



+Sn,k+1 [Ckd{Ck+l) + SkdiSk+l)] j{L0,L0), 
^{ifk - Vk+l) ■ 



COS 



1 1 

^ 2 (z-w)2 

H — 9 (: cos V2Tr{Lpk - (fk+i) :) (a;,w), 

Z — UJ \ L J / 

(A26) 

and similarly for the second term of Vsd, which proves finally 



Eq. (A25i. 



3. Free-field representation of tlie TIM and Z3 Potts CFTs 

In this section, we give a proof of the non-trivial identifica- 
tion pT] ) which occurs in triple-tube systems. 

It is first important to observe that the free massless bosonic 
field $2 has a very special compactification radius R2 — 
\/3Jtt in the classification of the CFT with central charge 
c = 1. At this radius, it displays a CFT with an extended 
symmetry: a TV = 2 (respectively iV = 1) superconfor- 
mal field theory (SCFT) when the bosonic field $2 is com- 
pactified along a circle (respectively an orbifold $2 ^ —$2) 
with radius i?2!^ Forgetting the contribution of the Majorana 
fermion the conformal symmetry of the non-interacting 
limit of model ( [56] l is Ising x [c ~ 1 SCFT], with central 
charge c — 3/2. The latter CFT can also be described in terms 
of the product of TIM x Z3 Potts CFTs with central charge 
c = 7/10 + 4/5 = 3/2. The precise conformal embedding 
has been derived by the authors of Ref|45] 

Z2 X (c = 1 iV = 1 SCFT) = P[MiX M5] ,(A27) 

where the projection P restricts the primary operators of 
M4XM5 to the subset: {^i*}<i>i^}}, (1 < r < p-l, 1 < 
s < p) being the primary operator of A4p CFT. 

As is well-known, the (left) stress-energy tensor of the Z2 
(or Ising) CFT expresses directly in terms of the Majorana 
fermion Tj = —tt : ^j^dS,}^ :. It is indeed straightforward 
to check that Tj satisfies the defining relation ( A8 1 of a stress- 
energy tensor with c— 1/2. The left Majorana is a primary 
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field with holomorphic weight h = 1/2: 
ei(0) , 9^(0) 



T/Wei (0) 
?i Wei (0) 



2z2 

1 



r/(0) + — aT/(0).(A28) 

ZTT 



For the = 1 SCFT with central charge c = 1, the stress- 



energy tensor is simply given by Eq. (A9i: Tq = — 27r : 
(9<i>2L)^ :. On top of this bosonic tensor, the iV = 1 SCFT is 
characterized by the existence of a fermionic current G with 
holomorphic weight h = 3/2*23 



To (z) G (0) 
G(z)G(0) 



3G (0) dG (0) 



2z2 



3To(0) , 39ro(0) 



(A29) 



In the particular case of the = 1 SCFT with central charge 
c = 1, the current G has a simple free-field representation in 
terms of the free-massless boson ^2L^ 



G = \/2 : cos(%/12^$2l) 



(A30) 



It is indeed easy to check, using Eq. \A5\ , that the OPEs ( A29 1 
are reproduced from this identification. 

From all these definitions, we find a new free-field rep- 
resentation of the left stress-energy tensors of the TIM and 
Z3 Potts CFTs in terms of a Majorana fermion and a free- 
massless boson: 



Tw 



\ti + ^To + : cos(\/12^$2l) (A31) 



where Ttim and T^g denote respectively the stress-energy 
tensor of the TIM and Z3 Potts CFTs. Using the results ( [A^ 
A29|, one can indeed show that Ttim (respectively T%^) satis- 



fies the definition ( A8 1 with c = 7/10 (respectively c = 4/5), 
together with the decoupling: Ttim (2)223 (0) ~ In addi- 
tion, we note that we have Ttim +7z3 = T/ +To in full agree- 
ment with the existence of the conformal embedding ( |A27[ ). 
We are now in position with the free-field representation 



(A31 1 to show the identification (57i. In this respect, let us 



introduce the following operators: 

0\ = fj-i : cos ^V37r$2^ : +0"! : cos 



(%/3^e2) : 



O2 = 



/xi : cos ^-\/37r$2^ : — ci : cos (^VSttQ2^ ■ ■(A32) 



Let us first show that Oi is a singlet under the Z3 Potts CFT 
and a primary operator of the TIM CFT with holomorphic 



weight h = 7/16. To this end, we need the following OPEs 
for a massless boson field that can be obtained from Eq. ( |A5| i: 



: cos(V127r$2L) : (z) : cos(V37r$2) : (0,0) 

g-i37r/4 



2z3/2 



cos(V37re2) : +2z : acos(V37re2) : (0,0) 



: cos(\/12^$2l) : (z) : cos(\/3^e2) : (0,0) - 

gi37r/4 |- -, 

[: cos(V3^$2) : +2z : 9cos(\/3^$2) :J (0,0) . 
^ (A33) 
The following Ising OPEs, which can be found in Ref|72j are 
also needed: 



{z)a, (0,0) 
(z)//i (0,0)^ 



2^7^21/2 

g-i7r/4 

2V^zi/2 



[/ii+4z9/ii] (0,0) 

[o-i +4zacri] (0,0). (A34) 



With all these identifications, it becomes straightforward to 
obtain the following result directly from Eq. ( |A31| i: 



Ttim (2)Oi (0,0) 
Tz3(z)Oi(0,0) 



701(0,0) 9Oi(0,0) 



16z2 



0, 



(A35) 



which states that Oi is indeed a primary operator of the TIM 
CFT with scaling dimension 7/8. It is thus proportional to the 
subleading magnetization CTtim- Using the correlation (50i, 



one can fix the normalization factor by evaluating the two- 
point function of Oi: 



(Oi(z,z)Oi(0,0)) 



(^TIm(2j ^) ^TIM 

(0,0)), (A36) 



so that we deduce Oi ^ cttim- Finally, from its definition 
( A32 1, it is easy to check that O2 is a primary operator of the 
Mi X M5 CFT with scaling dimension 7/8: 



(T 



TIM 



-rz3)(z)02(o,o) 



(T/ + To)(z)02(0,0) (A37) 
702(0,0) ^ 902(0,0) 
16z2 Z 



The P[M4 X M^] CFT contains only two primary operators 
with scaling dimension 7/8, i.e., drj^iy^ and CTiMezs- We thus 
find that O2 ^ crTiMfZa since Oi ^ (Ttim- 
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